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We have revisited the magnetic penetration depth data on the electron-doped cuprates 
P^-zCe^CuCU-y and La2- a: Ce a ;Cu04_, / . It is proposed that the transition between nodal-gap- 
like and nodeless-gap-like behavior upon electron-doping [e.g., see M. Kim et al., Phys. Rev. Lett. 
91, 87001 (2003)] can be due to a scattering of the quasiparticles in the d-wave superconducting 
state by an incipient or weak antiferromagnetic spin-density-wave. This conjecture is supported 
by the inelastic neutron scattering and angle- resolved photoemission experiments on some closely 
related electron-doped cuprates. 

PACS numbers: 74.72.-h, 74.20.Rp, 74.20.Mn 



I. INTRODUCTION 

A common opening sentence in papers related to the 
high-T c cuprates is, "After fifteen years of its discov- 
ery, the mechanism of superconductivity in the high-T c 
cuprates is still under debate...". Meanwhile, an inti- 
mately related question, the symmetry of the supercon- 
ductivity order also still possess some ambiguity, and 
nailing down the gap symmetry can lay a tread for at- 
tempts to unveil the superconducting (SC) mechanism. 
Today's situation is that, in the hole-doped (ft,-doped) 
compounds, tctragonally symmetric d a .2_ J) 2-wavc super- 
conductivity (dSC) is obtained unanimously in the un- 
dcrdoped and optimally-doped compounds^. In the over- 
doped compounds, diverged results like fully gapped 



d+is, d+id— or orthorhombic d r 



be ob- 



_ y 2+3ii may 

tained by different probes. For the electron-doped (e- 
doped) compounds, observations are more converged^. A 
nodal SC gap, most probably d x 2_ y £, is observed at un- 
dcrdoped. Upon electron-doping, a nodclcss gap is devel- 
oped as seen in the magnetic penetration depth (MPD) 
measurementSi^iSiiSiii. Moreover, the gap may become 
nodal again at overdopedii. The symmetry of the node- 
less gap is controversial and various possibilities have 
been proposecU2iiiii4, such as a conventional s-wave, or 
some unconventional symmetry waves like p+ip— . If 
one admits that the symmetry of the e-doped compounds 
could indeed differ from its /i-doped partner and trans- 
mutation of gap symmetry does occur, such a situation 
will be bewildering. It seems that there are different 
mechanisms of superconductivity governing the e-doped 
and /i-doped cuprates. Moreover, a quantum phase tran- 
sition (QPT) or more QPTs may be taking place in the 
e-doped cuprates. We think that the diverging observa- 
tions of the /i-doped cuprates at the overdoped regime are 
due to difficulties in the individual experiments, but the 
more converged result of the e-doped compounds needs 
to be seriously look into. 

More often the transition between nodal-gap and 
nodeless-gap behavior in the MPD, on which most con- 



clusions of a gap symmetry transition are based, is as- 
cribed to a QPT of some origini^. In this paper, we 
have chosen a more conservative but phenomenological 
approach. We take two observational facts as the input. 
A caveat beforehand is that those observations are not on 
the compounds which the MPD data were taken, due to 
the fact that a particular experimental probe usually does 
not access different compounds equally well. Neverthe- 
less, those observations are on the e-doped compounds. 
The first observation is due to the inelastic neutron scat- 
tering (INS) experiment, by which a soft-frequency and 
momentum q = Q = (tt, 7t) (commensurate) scattering is 
observed in the Nd2- a; Ce 2; Cu04 compound a 16 i 17 i 18 . The 
presence of such a soft excitation in contrast to the h- 
doped cuprates is in congruence with the fact that the an- 
tiferromagnetic (AF) phase is more robust in the e-doped 
cuprates, and destroyed only at higher doping levels^. 
The second observation is due to the angle-resolved pho- 
toemission spectroscopy (ARPES) on the Nd compounds. 
A pseudogap is seen to open up at the intersections of 
the AF magnetic Brillouin zone (MBZ) boundary and 
the Fermi surface (FS)S!l. The scattering that is respon- 
sible for the pseudogap thus has momentum (tt, tt) and is 
sharp in momentum space. Putting the INS and ARPES 
observations together, we propose that in those e-doped 
cuprates where the MPD data were taken, there also ex- 
ists a soft magnetic fluctuation with momentum (tt, tt) 
and it is coupled to the quasiparticles. 



II. FORMULATION AND RESULT 

We will study a model system with dSC co-existing 
with a weak AF spin-density wave (SDW). The weak 
SDW order is used to model an incipient SDW, or a 
soft-energy, momentum q = Q magnetic fluctuation. An 
implicit assumption in this approximation is, all ver- 
tex corrections in the current-current correlation func- 
tion (which is related to the MPD) are neglected, and 
only the self-energy which directly affects the quasipar- 
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tide (QP) spectrum is considered. As shown by Kampf 
and Schrieffe»2i, the change in the single-particle spec- 
trum due to a static SDW was essentially the same as 
that due to the scattering by a long coherence length, soft 
magnetic fluctuation (see also Ref. 122^) . Thus this mod- 
eling should be physically sound for the current problem. 

We start from a hamiltonian with dSC and AF SDW 
order, 



H = E & c L c k* + E[ AkC k T c -ki + H.c] 



E* 



(1) 



4i'cos k x cos k y — fx, 



where £k = — 2i(cos k x + cos k y 
t is the nearest-neighbor (NN) hopping which is set to 
1, t' = —0.4 is the next-nearest-neighbor (NNN) hop- 
ping used to reproduce a typical bowed FS, Q = (7r,7r), 
a = +1 ( — 1) for spin up (down), fj, is the chemical po- 
tential, and Ak and <& are the SC and SDW order re- 
spectively. For a particular band filling, a fixed a is used 
at all temperatures, superconducting and normal states. 
The slight change of y, is below one percent in all our 
cases, and is neglected. All equations in this paper are 
derived for a complex- valued Ak and a real- valued <£>, 



but we will be discussing only the d x i _ y 2 -symmetry SC 
state, i.e., Ak = A k = Agc(cos k x — cos k y )/2. The 
values of Age and <3? in our context are not solved in the 
usual mean field manner, but they are chosen. We have 
chosen Asc = 0.1 and $ = 0.3^ throughout this paper, 
and the resulting system is metallic (in accord with the 
real materials). 

Defining a Nambu operator 



Ck| 

c f 

C -kl 
C k+QT 

-k-Ql 



(2) 



one can then define a 4x4 Matsubara Green func- 
tion matrix G(k, t) component-wise by [G(k, T)] a p = 
— (T^kaC?")^/?)- The temporal-Fourier-transformed 
SC Green function matrix G (k, iuj n ) is found to be 



G b °(k,i w „) 



V 



<(k) 



r/.a 



4^ iu> n - E n<7 (k) ' 



(3) 



where af^(k) is the hermitian spectral weight matrix 
whose components are 



[< C (k)] n 
[^(k)] 

[«£00] 
k^(k)] 
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l.'i 



2\E + (kf - E^{kf] 

a 

2\E + (k)i - E^ky}E v<J (k) 
er$ f 1 



^a(k) 



[^(k) 



E va (k) 



(£k + 6c+q) f , 



{A k [S, CT (k) 2 - El +Q ] + A k+Q $ 2 }, 

[E v<y (k) 2 - $ 2 + £k£k+ Q + A k A k+Q i + (a + ek+ Q ; 
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2[£+(k) 2 -£_(k) 2 ] \E n „{k) 

2[ E+ {ky -L(k)2]^ ff( k) [ ^ (k)(Ak + Ak+Q) + €kAk+Q - Ck+QAk 



(4) 



r 



[a^(k)] 22 , [o^(k)] 2 3, and [a^(k)] 24 are obtained from [a^(k)] 44 are obtained from [a*°(k)] n , [a^(k)] 12 , 

[a^(k)]u, [a^(k)]i4, and {af£(k)} 13 respectively by and [a^(k)] 22 respectively by the exchange £ k <-> 

the substitution £ k — > -£k,£k+Q — > -£k+Q,A k — > £k+Q, A k <-» A k+ Q. The poles are at E va (k) = r)E a (k) 

A k ,A k+Q -► A* +Q ; while [a^(k)] 33 , [o^(k)] M , and (n, a = ±1), with 



^(k) 



\ 



^k+Q 



$2 



<J> 2 [(^ + ^ +Q ) 2 + |Ak + A k+Q |2], 



(•5) 



where i?k = (£ k + Ak] 2 ) 1 / 2 . Note that for the current vanishes, but is kept for better illustration. Each SC 
case Ak+Q = — Ak and the term |Ak + Ak+cj| in E a (k) branch is splitted by the SDW or vice versa speaking, 
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FIG. 1: Normal state DOS with SDW order $ = (black), 
0.2 (blue), and 1 (red) are shown in (a). The chemical po- 
tentials are set at zero. We will use $ = 0.2 for all calcula- 
tions later in this paper. The near-Fermi-level spectral weight 
A(k) (for <E> = 0.2 and Asc = 0.1) are shown for the cases (b) 
fj, = -0.3 and (c) +0.1. 



forming four branches of excitations E va (k). 

A few more words on what we mean by a weak SDW 
is desired. Figure da) shows the normal state DOS for 
cases of different magnitudes of the SDW order, and it 
is seen that for our band with a weak SDW, no gap is 
formed at any chemical potential in contrast to the strong 
SDW case. Nevertheless, the FS may still be fragmented 
at some band fillings as seen in the distribution of near- 
Fermi-level spectral weight [Fig. IHb) and (c)] . The near- 
Fermi-level spectral weight as a function of momentum is 
defined by A(k) = f^^j du; Tr^ImGff (k, w+i0)f(w), 
where f(ui) is the Fermi distribution function but here it 
is taken as the Heaviside step function 9(— u>) for T = 
limit. The integration limit "0.15|t|" is an arbitrarily cho- 
sen small parameter, which is larger than Age however. 

We now examine the low-temperature behavior of 
MPD. The MPD tensor X^T) is related to the zero- 
momentum and zero-frequency current-current correla- 
tion function. The current operator can be written as 



U(q) = -e *k7„(k+!)*k +q , (6) 

kGMBZ 
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FIG. 2: Low-temperature magnetic penetration depth at var- 
ious band fillings is plotted as A(0) 2 /A(T) 2 versus the reduced 
temperature T/T c (2 Asc /T c = 4 is used). The chemical po- 
tentials are fi = —0.3, —0.1, +0.1, +0.3 (fillings from 1.16 to 
1.30 for our band). 



where 7 M (k) = Diag[v£l 2 , v£ +Q l 2 ], v 11 = <9£ k /<9k M , 1 2 = 
Diag[l,l] is the 2x2 unit matrix, e is the electronic 
charge, and k is sum over the AF MBZ. A^(T) is then 
given by 

1 



E 



Tr [G sc (k, uu n ) lfi (k)G sc (k, iw n )7„ (k) 

- G ,N (k, ^, 1 )7M(k)G N (k,i w „) 7 ,(k)] , (7) 

where G N is the normal-state Green function matrix, 
which can be simply obtained from G sc by setting Age 
to zero, but keeping the same SDW order <3>. Summing 
off the Matsubara frequencies, Eq. (J7J is reduced to 



1 



/[£V(k)]-/[iWk)] 



\^{Tf N E l2 £^(k)- 

^ k 7i<T7i' cr ' 

xTr [a^(k) 7 „(k)o^(k)7„(k) - (A sc 



-0)]. (8) 



Due to the tetragonal symmetry, \~„ = 5^ v XZz, and 
X xx = X yy = A. Some further simplification can be made 
and eventually this integral is evaluated numerically. The 
main contribution to the integral comes from the large 
[f{ijj) — f(u>')]/[u) — a/], i.e., near-zero frequency regime, 
and it is therefore a sum over low-energy thermal excita- 
tions. When lu = lu' (e.g., when r\ = rf and a = a'), the 
quantity [f(u>) — f(u)')]/[u) — u>'] is defined as the limit 
df /du. 

Figure |21 shows the low-temperature MPD. The same 
Age and $ are used within the temperature range con- 
sidered, and this approximation should be valid at least 



4 



0.008 



0.006 



0.004 



O 



0.002 



0.000 





1 1 


i i i i 










\l= -0.3 




- 


|i= -0.1 


■ ^s" — 




H = + 0.1 


s^. . ' 




H = + 0.3 


— . ' 










s - ' 




> 






j/s' . ■ 
K*' , - * 


- 









0.00 0.01 0.02 0.03 0.04 0.05 0.06 

T 2 /T 2 




FIG. 3: The specific heat CV is plotted versus T 2 . The band 
parameters are the same as those in Fig. [2] 

at the very low temperature regime, say below ~ T c /A. 
At approaching zero temperature, we find linearly de- 
creasing behavior A(0) -2 — A(T)~ 2 oc T at low and high 
band fillings, which is essentially due to a nodal dSC; and 
negligible temperature dependence at some intermediate 
band fillings, which is due to a nodeless excitation gap in 
the single particle spectrum (see Fig.EJ. 

The transition between nodal-gap and full-gap be- 
havior in the MPD is compared to the results of low- 
temperature specific heat. Figure [3] shows the constant 
volume specific heat 

C ^ T ) = E E na (k)f[E va (k)} 

keMBZ.ijff 

= 2 E ^/[^(k)]/[-^(k)], (9) 

keMBZ,er 

where df(e)/dT = (e/T 2 )/(e)/(— e) is used in the second 
line. In congruence with the MPD, the T 2 dependence 
of C v (T) at low temperatures can be ascribed to the lin- 
ear energy-dependence of the low-energy density of states 
(DOS) of the dSC (i.e., N(uj) oc |w| at u ~ 0, see Fig.gJ; 
and the case of negligible tcmpcraturc-dcpcndcncc is as- 
cribablc to a nodeless gap (see Fig.0J. 

Figure H shows the DOS [i.e., N(uj) = 
-(7r^)- 1 J] k Im[G sc (k,cj + i0)] n at different band 
fillings. The full gap24 that appears at intermediate 
band fillings explains the small temperature dependence 
of A (see Fig- EJ an d C v (see Fig. [2Jl at low temperatures. 
Interestingly, a subgap twin-peak structure 2 ^ is seen in 
the case of the MBZ boundary cuts the FS but the gap 
remain nodal-like (the case of fj, = —0.3 in Fig.0J. Some 
"peak-dip- hump" structure may also appear (e.g., see 
the negative energy side of the fj, = —0.1 case in Fig.0J. 

The formation of the full gap when the FS is close 



FIG. 4: The DOS in the SC state near the Fermi level is 
shown with chemical potentials /j, — —0.3, —0.1, +0.1, and 
+0.3. The inset shows the Fermi surfaces when there was no 
SDW order, at chemical potentials n — —0.3, 0, +0.3. The 
MBZ boundary is also shown. 

to the MBZ boundary is illucidated in Fig. |SJ Figure [S] 
shows the single- particle spectrum — 7r _1 Im[G sc (k, u + 
i0)]n along the diagonal-k directions, where a plain 
d a .2_ y 2-wave SC system would have its gapless single- 
particle excitations. At the presence of a SDW, the 
diagonal-k QP peaks do not disperse to the Fermi level 
when the FS (defined at the absence of SDW) gets close 
to the MBZ boundary [Fig.^b). See also Fig.^c)], and 
this explains the full gaps in Fig. 0] Scattering by the 
SDW has opened up a gap at the vicinity of the MBZ 
boundary. Interestingly, since the ARPES can observe 
only the negative energy part of the single-particle spec- 
trum, the QPs at the SC nodes that are gapped by a SDW 
[see Fig. Elb)] may be seen to disperse like QPs that are 
gapped by a SC gap. Moreover, though the DOS may 
not vanish on the Fermi level, a forbidened window in 
w-space at below the Fermi level may still be found for 
the diagonal-k QPs [see Figs. 0] and EJc)]. 



III. CONCLUDING REMARKS 

We have used the same Asd T c , and $ throughout 
our discussion for simplicity. Taking into account the 
differences of those parameters in various cases should 
make no qualitative difference in our conclusions as we 
are essentially discussing thermally activated behaviors 
due to a full gap, or scale-less behaviors due to a nodal- 
like gap. As long as <!> is finite, there exists a range of 
doping (within the realistic doping regime) that can give 
rise to a full gap. Furthermore, we also believe that our 
conclusion will remain qualitatively the same when the 
vertex correction is taken into account, since the gap at 
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FIG. 5: Single-particle spectrum -7r _1 Im [G s (k,w + i0)]n 
along the diagonal-k direction is plotted from (0.457T, 0.457r) 
(lower most) to (0.55-7T, 0.557r) (upper most), for (a) fi — —0.3, 
(b) /Li = —0.1, and (c) = +0.3. The heights of the spikes 
are drawn proportional to the weight of the poles. 

the MBZ boundary will always be opened up due to the 
mixing of QPs at k and k + Q. As the would-be-gaplcss 
points on the FS get close to the MBZ boundary, they 
are pushed to higher energies and a full gap in the total 
DOS is formed. 

We have adopted a phenomenological approach rather 
than obtaining the SDW order from a self-consistent 
mean- field-theory (MFT). This is because naive MFTs 
are usually insufficient to describe the very weak (prob- 
ably incipient) but persisting SDW indicated by the ex- 
periments. In such a situation, the location of the FS but 
not the size of the SDW, is the dominant factor in deter- 
mining the low-energy excitations in the dSC state. Note 
that though we have assumed a static SDW, the result is 
in fact equally applicable to the case of an incipient SDW 
or a soft AF spin fluctuation, as the electrons around the 



Fermi level arc "fast" modes, and any relatively "slow" 
spin fluctuations are essentially static. 

Though the properties of the dSC state can be dras- 
tically modified at some of the band fillings, the normal 
state properties are expected to remain normally. For in- 
stance, the resistivity will always have a quadratic depen- 
dence on temperature because in the normal state there 
exists no gap in the DOS near the Fermi level [Fig.^a)]. 
On the impurity scattering effect on the MPD, the linear 
temperature dependence of A(T)~ 2 at low temperatures 
in the clean limit that occurs when the d-wave nodes on 
the FS are away from the MBZ boundary, is expected to 
cross over to quadratic in T (in the resonant scattering 
limit regardless of the existence of the SDW or not. 
The reason for the independence of SDW is that in those 
cases the low-energy nodal QPs are essentially unaffected 
— a fact that can be seen in the single-particle spectral 
function [see Figs. EJa) and (c)], and manifested in the 
linear in \w\ dependence of the DOS at low energies (see 
Fig. 0J. As the low-temperature MPD is related only 
to the low-energy excitations of those unaffected QPs, 
a theory including the effect of impurity scattering will 
be essentially the same as that without SDW, which is 
already widely available in literature 26 . 

At the overdoped compounds, two possible factors 
might account for the recovering of a nodal gap as seen 
in Ref. The first is, the diagonal portions of the FS 
have proceeded to the other side of, and leave the MBZ 
boundary (as studied in the current paper). The second 
is, the strength of the magnetic fluctuation is weakened. 
Precise ARPES and INS measurement on the compounds 
may be used to resolve the possibilities. 

We have investigated the transmutation between 
nodal-gap- like and nodeless-gap-like behaviors in some of 
the electron-doped cuprates. A scenario inspired by ex- 
periments on some related compounds is proposed. We 
propose that there exists a weak or incipient SDW order 
in the compounds, and it is coupled to the QPs. This 
scattering can push the would-be-gapless diagonal QPs 
to higher energies, and a d K 2_^2 -symmetry SC system 
can behave like it is nodeless. 
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